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Abstract. We study the conductance of mesoscopic graphene rings in the presence 
of a perpendicular magnetic field by means of numerical calculations based on a tight- 
binding model. First, we consider the magnetoconductance of such rings and observe 
the Aharonov-Bohm effect. We investigate different regimes of the magnetic flux up 
to the quantum Hall regime, where the Aharonov-Bohm oscillations are suppressed. 
Results for both clean (ballistic) and disordered (diffusive) rings are presented. Second, 
we study rings with smooth mass boundary that are weakly coupled to leads. We show 
that the valley degeneracy of the eigenstates in closed graphene rings can be lifted by 
a small magnetic flux, and that this lifting can be observed in the transport properties 
of the system. 



1. Introduction 

Since their experimental discovery, graphite monolayers, also known as graphene, have 
attracted a huge amount of interest among both experimentalists and theorists due to 
the linear low energy dispersion and various properties stemming from this unusual 
dispersion [HE]- For instance, graphene has opened new perspectives for mesoscopic 
physics, such as pseudo diffusive conductance at the Dirac point [31 HI [5], specular 
Andreev reflection [6], or the signatures of symmetries of the graphene Hamiltonian 
in the conductance of diffusive wires [3 El El [iDl [H] and ballistic cavities [121 113] • 

Mesoscopic rings can be considered as prototypical devices in mesoscopic physics, 
as they show one of the most basic coherence effects, namely the Aharonov-Bohm 
(AB) effect [T^ IT5l [T6]: oscillations of the transmission, or dimensionless conductance, 
T = {h/2e^)G as a function of the magnetic flux $ through the ring. The reason 
for these oscillations is the phase difference Aip = 27t^/^q between electrons travelling 
along the different arms of the ring. Here, = h/e is the magnetic flux quantum. Fifty 
years after its discovery, the AB effect is one of the most well-known manifestations of 
quantum interference within and beyond the field of condensed matter physics. 

Hence it is rather surprising that the AB effect in graphene has up to now received 
only a little attention in the literature. Most notably, there are only two experiments 
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on graphene AB rings manufactured by electron beam lithography [T71 [TS] , one of them 
leaving many open questions on the physical origin of some of the observed effects |17j . 
From the theory side, there is only one numerical study of the AB effect in graphene 
rings; it focuses on the effects of valley-polarized currents, i.e. on the few-mode or low- 
doping regime in the leads [19j. In this work, we will in contrast also consider the 
many-mode or high-doping regime. 

In addition to these studies focusing on the transport properties of open rings, there 
has been a proposal to use the Aharonov-Bohm effect in closed rings to form qubits: The 
energy spectrum of a closed graphene ring with infinite mass boundary condition [20] 
has been calculated in reference [21] , where the authors find that the valley degeneracy 
of the energy levels is lifted as soon as a magnetic flux pierces the ring. This effect has 
also been found for chaotic graphene quantum dots [13]. Note that this aspect is not 
present in AB rings realized e.g. in semiconductor heterostructures and metals. It is 
connected to a special antiunitary symmetry of the Dirac Hamiltonian, which describes 
graphene well for low Fermi energies. In this work, we will show that the lifting of the 
valley degeneracy is also visible in the transport properties of graphene rings. 

The paper is organized as follows: In the first part, we investigate the AB effect 
of graphene structures by numerically calculating the transmission of rings attached to 
infinitely extended leads. We study both small rings in the one-mode regime and large 
rings with many modes propagating in both the leads and arms of the ring. In the latter 
we especially consider the high-field regime and the effects of disorder. In the second 
part of this work, we show that the breaking of valley- degeneracy by a magnetic field 
is also visible in the transport properties of graphene rings. We do this by numerically 
calculating the transmission of graphene rings that are weakly coupled to two leads. 
This transmission shows peaks as a function of the Fermi energy Ep which correspond 
to the energy levels of a closed ring; the lifting of their degeneracy can be observed as a 
splitting of the transmission peaks upon applying a magnetic field perpendicular to the 
ring. 

For our numerical work, we use a tight binding model taking into account the 
2p2-orbitals of the carbon atoms, leading to the Hamiltonian 



with i and j beeing nearest neighbor sites in the first sum. The magnetic field is included 
via the Peierls substitution t^- = — texp [i^ f^' A - dr). The second term accounts 



for a staggered on-site potential, i.e. mj = m{ri) is positive (negative) if is located 
on sublattice A (B). Such a staggered potential corresponds to a mass term in the 
effective Dirac equation and will be used in the second part of this paper to suppress 
the inter- valley scattering that breaks the valley degeneracy [13]. The lattice points 
are determined by cutting a ring out of the graphene sheet [cf. figure [11(b)] . In order 
to solve the transport problem to obtain the dimensionless conductance T within this 
tight-binding model, we use an adaptive recursive Green function method [22] . 




(1) 





Figure 1. (a) Schematic of the systems studied numerically. The parameters defining 
the shape are the inner and outer radius i?i and R2, respectively, and the width Wl of 
the infinitely extended leads. The dashed line marks the points where the mass term 
used in Section [3] is zero, (b) Example of a system defined by cutting a ring out of a 
graphene sheet [the system was used for figure [5] (a) and (c)]. 



2. Aharonov-Bohm effect in graphene rings 

We now investigate tlie transmission of grapfiene rings witli two attaclied leads under 
tlie influence of a perpendicular magnetic field B, giving rise to a flux $ through the 
ring. In the following we define ^ = B ttR"^ as the flux through a circle with the average 
ring radius R = + -R2)/2 (cf. figured]). Because electrons can travel along one or 
the other arm of the ring, the transmission T($) is expected to oscillate with periods 
A$ = ^Q/n,n G N, as mentioned in the introduction. The reason why more than one 
single period, namely $0) may occur is that the electrons do not necessarily have to 
leave the ring after traversing an arm once, but rather they may circle around the ring 
several times before exiting, giving rise to higher harmonics with n > 1. 

We begin by considering small graphene rings in the one-mode regime. In figure 
[2](a) and (b) we show the numerically calculated transmission as a function of $ for 
two small rings with radii of about 10 nm and widths of about 3nm (see figure caption 
for details). One of the rings has zigzag type leads while the other has armchair type 
leads. For both, one can clearly see the AB oscillations with period $ = $o- To expose 
a few more details of the frequency content of our data, we show the power spectra 
of the AB oscillations in figure [2](c) and (d). There we find pronounced peaks for the 
fundamental frequency and for several higher harmonics (we show only the first five). 
Note that we plot the spectra on a logarithmic scale, since the fundamental frequency 
is strongly dominating over the higher ones. The peaks in the spectrum lie very close 
to multiple values of l/$o- Because of the finite width of the rings, one has certain 
allowed deviations from these values. For example, for the system of figure [2](a), the 
fundamental frequency is expected to lie between 0.74/$o and l.S/^o, obtained from 
the inner and outer ring radii Ri and R2, respectively. The fact that the peaks are 
slightly shifted to frequencies lower than n/$o means that the effective radius of our 
rings is slightly smaller than the average radius R. We also performed calculations for 
various systems with different parameters and found the same behavior irrespective of 
the details like radii, width or lead type. 
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Figure 2. Aharonov-Bohm oscillations [(a),(b)] and corresponding power spectra 
[(c), (d)] of small rings in the one-mode regime. (a),(c) Ring with zigzag type leads and 
Ri — 9.84 nm, i?2 = 13.0 nm, Wl = 3.20 nm. (b), (d) Ring with armchair type leads 
and Ri = 9.84 nm, i?2 = 12.5 nm, Wl = 4.92 nm. f{l/A^) is the Fourier transform 
of T(<I>) where A$ is the period of the oscillation related to the Fourier transform 
variable by fc = 27r/A$. 



Next, we turn to the case of larger graphene rings in the many-mode regime, i.e. for 
high values of doping. This is the regime applicable to the available experiments of 
Refs. p!7| [T8] . In the numerical simulations we consider rings with an average radius 
R 55 nm. While this size is still smaller than in the experiments {R ~ 500 nm in 
Ref. 1^7\ and R ~ 250-300 nm in Ref. [IH]), the rings contain more than a hunded 
thousand atoms and hence can be considered as mesoscopic objects. Therefore, we do 
not expect that larger rings would show any fundamentally different physical properties 
than the systems considered here. 

In figure [3], we show the magnetoconductance of a large ring for a wide range of the 
magnetic flux, up to the quantum Hall regime. Up to a flux of about 300 $0) we see AB 
oscillations with amplitudes approximately in the range AT ^ 0.5 — 1 [see also figure 
|5](a)]. For fluxes larger than 300 $0; the oscillation amplitude becomes significantly 
smaller (AT < 0.5), while for even higher values of $ the oscillations vanish completely 
and the transmission is equal to the maximum value given by the number of quantum 
Hall edge channels in the leads. This behavior can be understood by the following 
picture. The cyclotron diameter dc in (bulk) graphene is given by 

dc = —^, (2) 

where Ep is the Fermi energy and vp is the Fermi velocity. In our case Ep = 0.3 1 and 
the width of the ring arms is ~ 18 nm. This means, that at a flux of $ ~ 230 $0) 
the cyclotron diameter is equal to an arm's width. When $ becomes comparable to this 
value, it becomes more probable that the electron enters one arm than the other, because 
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Figure 3. Magnetotransmission of a large ring with i?i = 45.5 nm, i?2 = 63.7 nm 
and Wl — 13.5 nm for a wide range of the magnetic flux $ (sohd black line) in the 
multi-mode regime. The red dashed line shows the number of occupied transverse 
channels in the leads. The insets show power spectra for different regimes: left inset 
$ = 0-45 $0, right inset $ = 305-350 Mo- 



tile magnetic field dictates the direction of classical propagation. Therefore, interference 
is reduced. When the magnetic field gets larger and larger, it becomes more and more 
unlikely for the electron to enter the second arm, and finally the interference vanishes 
when the quantum Hall edge channels fully form. Since backscattering is also strongly 
reduced at large magnetic fields, the transmission is only limited by the number of 
propagating lead channels. This number however is reduced for a fixed Fermi energy 
as the magnetic field is increased (the distance between the Landau levels increases 
with B), which leads to the step-like transmission in this regime. Furthermore, classical 
cyclotron effects are also important at the entrance to the ring. At low fields, a large 
part of the wavefunction entering the ring is simply reflected back into the lead by the 
inner wall of the ring. As the field gets higher, the probability to turn around the corner 
into one of the arms increases. 

This scenario is shown in figure HI the absolute value of the current density in 
the ring is given for values of the flux $ corresponding to the different regimes. For 
$ = 113 $0 [flgure |4](a)] - the low transmission regime - the current flows more or 
less equally through both arms. As the flux is increased, the electrons are forced more 
and more into the upper arm, and at very high fluxes, when electrons are essentially 
travelling only along the upper arm, interference gets suppressed and flnally vanishes. 
Note that also the quantum Hall edge channels become clearly visible at high fluxes. It 
is worth noting that the small AB oscillations in the regime where $ > 300 $0 are not 
due to electrons entering the lower arm of the ring from the left lead. If this were true. 
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Figure 4. Absolute value of the current density in the ring of figure [31 The electrons 
enter the ring from the left lead. Blue means low and red means high current density, 
(a) $ = 1 13 $0 (b) $ = 204 $o (c) * = 328 $o (d) $ = 453 $o- 

the current in figure |l](c) would be localized at the inner border of the lower arm, which 
is obviously not the case. In this regime, the interference is due to paths going along 
the upper arm and leaving the ring after a half circle interfering with paths which also 
start along the upper arm but which go around the ring one more time before exiting. 
Such paths also give rise to oscillations with period $o- Note that this effect is present 
only because the lead width in the chosen geometry is somewhat smaller than the width 
of the ring, as then the quantum Hall edge channels in the ring arms may be scattered 
at the lead opening and enter the opposite arm. Hence, this regime corresponds to 
interferometry with quantum Hall edge channels. When the lead width is larger than 
the ring arm width, scattering at the lead openings is suppressed, and no interference 
is observed. The power spectra shown as insets in figure [3] also support this picture: 
For low fields, several higher harmonics are present, whereas at high fields only the 
fundamental frequency is visible because multiple windings around the ring are strongly 
suppressed in the quantum Hall regime. 

The magnitude of the conductance seen in the quantum Hall regime - an odd integer 
times 2e^ /h (see figure [3]) - is consistent with the bulk quantum Hall effect [231 121] as 
well as with theoretical results obtained for graphene nanoribbons p^. (Recall that in 
a two terminal geometry the conductance shows the quantization of a^y) The steps of 
magnitude 2 are easily understood in terms of the valley degeneracy of the subbands in 
the leads. The quantization at odd rather than even integers is explained by a Berry 
phase effect involving the pseudospin associated with the sub-lattice degree of freedom 

PI- 

In the experimental work of Ref. a significant increase of the AB oscillation 
amplitude was observed at high fields when Wx ^ dc- Note that our data in figure [3] 
does not show any such increase (if anything, it shows the opposite). To investigate 
this further, figure [5] shows the oscillating part of the magnetoconductance for a range 
of $ near the value where the cyclotron diameter becomes equal to the width of a 
ring arm [around 230 $o according to ([2])]. The authors of [TTj speculate that an 
asymmetry between the two ring arms could be responsible for the increase of the 
oscillation amplitude. Our numerical calculations, however, do not confirm this: for 
figure [5](b) we reduced the width of one arm by 10 percent, but also in this case, an 
inrease of the oscillation amplitude at high magnetic fields is absent. 
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Figure 5. Oscillating part of the magnetotransmission. The average transmission of 
one period around the corresponding value of $ has been subtracted from T($) to 
obtain Tosc('I?)- Top panel: System of figure [3l Lower panel: Same system as in the 
top panel, but the outer radius of one ring arm is reduced from 63.7nm to 61.9 nm 
which results in a difference in the arm widths of 10 percent. 



The formation of quantum Hall edge channels in a clean system happens for 
dc ~ Wa, as observed in our simulations (figure [3]). In the data of Ref. [I7|, however, 
AB oscillations are still visible even for much higher magnetic fields, when dc < Wa- 
Furthermore, despite being ten times larger than the rings considered in our numerical 
studies, the average conductance in the highly doped regime is only on the order of a 

2 

few ^. These two features are a strong indication that these experiments are in the 
strongly disordered regime. 

Hence, we now turn to briefly discussing the effects of disorder on the AB 
oscillations. For this, we consider both edge disorder (employing the model of edge 
disorder from Ref. [26]) and bulk disorder (using Gaussian disorder as in Ref. |28j). 
Figure El^a) and (b) show the magnetoconductance for edge and bulk disorder, 
respectively. In the case of edge disorder only, we still observe an increase of the 
transmission with magnetic field due to classical cyclotron effects as in the clean case. In 
addition, the onset of quantized Hall conductance is only shifted slightly to higher fields. 
In fact, since the interior of the ring arms is ballistic, we still expect quantum Hall edge 
channels to form for dc ~ Wa- Moreover, as in the clean case we observe interference of 
quantum Hall edge states for $ > 300 $o- In contrast, for bulk disorder the transmission 
probability remains - apart from conductance fluctuations - virtually unchanged as the 
magnetic field increases. In this case, the disorder broadens the Landau levels enough 
to prevent the formation of edge channels. The similarly of this behavior to that of 
the experiments reported in Ref. [IT] suggests that the experiments are in the "dirty 
limit" , dominated by bulk disorder. Such a conclusion is also in agreement with Ref. [18] 
where their experimental results for AB rings are explained in terms of diffusive metallic 
transport. 

Further characteristics of the diffusive limit include that the system remains 
isotropic for higher fields, both ring arms carry an equal amount of current, and higher 
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Figure 6. Magnetotransmission of disordered rings with i?i — 45.5 nm, R2 = 63.7 nm 
and Wl — 13.5 nm in the multi-mode regime (solid black lines). The insets show 
corresponding power spectra, (a) Ring with edge disorder. Here, we employed the 
etching disorder model of Ref. [53] with A'swocps = 30 and pi = 0.05 for i = 1 ... 30 (for 
the details on the parameters of this model, see Ref. (b) Ring with bulk disorder. 

For the bulk disorder, we choose a Gaussian disorder potential with Ni„^p/Ntot = 0.05, 
6 = O.lt and f = 2a (for the details on the parameters of this model, see Ref. [27]). 
Note that although the Gaussian disorder potential is smooth on the lattice scale and 
hence does not scatter the valleys, valley scattering nevertheless happens at the abrupt 
termination of the lattice at the edges. 
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Figure 7. Oscillating part of the magnetotransmission for disordered rings. Top 
panel: System with edge disorder. Lower panel: System with Gaussian bulk disorder. 
Disorder parameters are as in figure [6l 



harmonics are suppressed in comparison to the clean case (see figure [3]). Finally, figure 
[7] shows the oscillating part of the magnetoconductance in the vincinity of ~ Wa- 
Again, the size of the oscillations remains essentially unchanged with increasing magnetic 
fiux. 

A closer inspection of the oscillations reveals that in the case of edge disorder, there 
are also aperiodic oscillations on top of the periodic AB oscillations for $ ^ 250 — 300 $0? 
where we observe the onset of edge channel formation. These aperiodic oscillations 
are presumably due to localization and resonant tunneling at the disordered edge, 
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since there is still backscattering among the not fully formed edge states. This 
interpretation is supported by the fact that the aperiodic oscillations are absent in 
the case of bulk disorder, where the edge channel formation is completely suppressed. 
Furthermore, note that aperiodic oscillations are also present (although to a less degree) 
for $ ^ 190 — 250 $0 in the clean rings, as the abrupt lattice termination also serves as 
a source of local disorder. 

3. Breaking the valley degeneracy in graphene rings 

In the following, we demonstrate numerically that the lifting of the valley degeneracy 
of eigenstates in rings by a magnetic field perpendicular to the ring is also visible 
in transport properties of the ring. In Ref. |21], the authors calculate the energy 
spectrum of a closed graphene ring with infinite mass boundary condition [20] within 
the effective Dirac theory. The energy levels E!^I^ are labeled by a radial quantum 
number n, an angular momentum quantum number m, valley index r and the sign of 
the energy. At zero magnetic field, the authors find from their calculation the degeneracy 
^n(-m) _ ^nm^ which is uot prcscut whcu a perpendicular magnetic field is applied. It is 
worth to point out that this is due to the fact that in addition to the usual time reversal 
(TR) operator T = {ay ® Ty)C, there is another antiunitary operator commuting with 
the Hamiltonian at zero magnetic field, namely % = {ay ® t^)C In contrast to T, % 
is associated with a symplectic symmetry, since = — 1. This ensures the degeneracy 
of both valleys (Kramer's degeneracy [29]). Note that this reasoning holds only if the 
boundary condition preserves the symmetry connected with 7^; such is the case for the 
infinite mass boundary condition used in Ref. [21] but is not the case, for example, 
for valley mixing armchair boundaries. Once a magnetic field pierces the ring, the 
Tv symmetry is broken and the valley degeneracy is lifted. This mechanism has been 
numerically observed in ballistic graphene quantum dots with smooth mass confinement 
[T3] , where the amplitude of universal conductance fluctuations was considered. Here we 
follow a different idea to probe the effect of breaking valley degeneracy with a magnetic 
field, namely we calculate numerically the conductance of a ring weakly coupled to two 
leads. The conductance is then expected to be peaked at values of the Fermi energy 
that match the energies of eigenstates of the closed ring. 

In order to have valley degenerate states at zero magnetic field, we have to avoid 
scattering at local armchair boundaries, as mentioned above. To model the infinite mass 
boundary numerically, we use smooth mass confinement as introduced in [13]. Due to 
the finite escape time of electrons injected into the ring from one lead, the effects of inter- 
valley scattering will be suppressed in such a system, provided the mass confinement is 
smooth enough. We use a mass term that is zero in the middle of the ring and increases 
quadratically towards the boundaries [cf. figure [T](a)]: m{x,y) = u"^ 6^{x,y)/2, where 
6{x, y) = a/x^ + — ^ for a ring with average radius R centered at x = y = 0. In the 
leads we use a consistent mass term, with S{x,y) = y for leads in the x-direction. In 
the region where the leads are attached, we match both terms in a smooth way. In this 
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Figure 8. Transmission of rings with smooth mass term {lu = 0.05\/i/a) and 
i?i — 9.84 nm, i?2 = 19.7 nm, Wl = 9.84 nm, with leads that are strongly (black) 
and weakly (blue) coupled to the system as a function of the Fermi energy in the 
two- mode regime (transverse lead channels are twofold valley degenerate). 



region we also introduce a mass- induced barrier, that can be varied in height, in order 
to tune the couphng strength between leads and ring. 

Figure [8] shows the transmission of rings with both strong and and weak coupling 
as a function of the Fermi energy Ep- We checked that the transmission is insensitive 
to the orientation of the system with respect to the lattice orientation at the energies 
considered. This means that the detailed edge structure of ring and leads do not play a 
role, and implies that the mass confinement is smooth enough to suppress inter-valley 
scattering. Upon reducing the coupling strength, the oscillations in the transmission of 
the ring with strongly coupled leads turn into sharp peaks. The system studied is a ring 
with Ri = 9.84 nm, i?2 = 19.7 nm and width Wl = 9.84 nm. Note that the effective 
width of leads and ring is reduced due to the mass confinement. The mass potential is 
0.5 1 at the physical borders of the system for the value of u used (see caption of figure 

ED. 

To find out, wether the energy spectrum of our ring shows the valley degeneracy 
connected with the symmetry described by 7^, we try to break this supposed degeneracy 
by applying a magnetic field perpendicular to the ring. Figure [9](a) shows the numerical 
results, where the magnetic flux through the area irR^ with R = 14.8 nm is $ = 
(dashed lines) and $ = 0.17 $0 (solid lines), respectively. We find that the transmission 
peaks split and that their heights are in general reduced. This is the behavior expected 
from degeneracy lifting, as the weight of each split peak now corresponds to only a 
single, non-degenerate level. Note that for E < O.lAt the peaks partially overlap; the 
splitting of the valley degeneracy is particularly clear for E > 0.14 1 where the peaks 
are well seperated, as shown in figure [9](b). These results show unambiguously that the 
lifting of the valley- degeneracy by a magnetic flux in weakly coupled graphene rings can 
also be observed in the transport properties. 




Figure 9. Transmission of the weakly coupled ring from figure [5] for zero flux (blue 
dashed lines) and a flux $ = 0.17 (red solid lines), (b) shows a zoom on the 6 (or 
12, respectively) rightmost peaks of (a). 



4. Conclusion 

In this work we liave investigated two different aspects of mesoscopic graphene rings: 
First, we have numerically studied the Aharonov-Bohm effect in rings cut out of the 
graphene lattice. Both small rings with a single mode as well as large rings supporting 
many modes exhibit clear h/e Aharonov-Bohm oscillations. Although disorder has a 
strong influence on the average transmission of graphene rings especially in the high 
field regime, the Aharonov-Bohm oscillations are influenced only very little by both 
edge and bulk disorder. Second, we have shown that the signature of the splitting of 
the valley- degeneracy by a magnetic field in rings with a mass confinement can also be 
observed in the transport properties of rings that are weakly coupled to leads. 
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